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OPTIMUM  PREDICTION  FOR  A 
TORQUE  SATURATING  SERVO 

SUMMARY 

This  study  was  made  to  examine  the  use  of  a  particular  type  of  nonlinear 
control  for  a  torque  saturating  second  order  servomechanism.    The  control 
involves  mechanization  with  diode  limiters  and  operational  amplifiers  of  a 
torque  switching  curve  in  the  phase  plane  of  error  and  error  rate,  with  a 
region  of  linear  operation  about  the  curve.    The  switching  curve  is  shaped  to 
give  prediction  for  optimum  time  response  and  essentially  no  overshoot  for 
step  and  ramp  inputs. 

The  results  of  simulation  on  an  electronic  differential  analyzer  show  a 
large  reduction  in  both  response  time  and  overshoot  for  a  servo  with  prediction 
as  compared  to  the  same  servo  with  proportional  plus  error  rate  control. 

The  introduction  of  an  approximation  to  the  true  switching  curve  and 
moderate  amounts  of  friction  had  little  influence  on  the  operation  of  the  pre- 
diction control.    The  shape  of  the  curve  could  be  altered  to  optimize  for  any 
particular  combination  of  Coulomb  or  viscous  friction  coefficients;  however,  the 
improvement  was  small. 

The  application  of  the  nonlinear  control  to  an  instrument  servo  consisting 
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of  a  400  cycle  per  second,  two-phase  motor  coupled  directly  to  a  one-turn  film 
potentiometer  substantiated  the  results  obtained  by  simulation.    The  character- 
istics of  the  signal  from  the  film  potentiometer  permitted  use  of  an  approximate 
differentiator  in  the  control  circuit.    With  this  system  a  resolution  of  one  part 
in  40,  000  and  a  response  time  of  0. 06  seconds  to  a  half- turn  step  input  were 
obtained. 
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SYMBOLS 

C  Coulomb  friction  torque,  lb.  -ft. 

Ce  error  rate  coefficient,  2  J/<s> 

e  control  voltage 

f  viscous  friction  coefficient,  lb.  -ft.  -sec. 

i  subscript  meaning  input 

I  moment  of  inertia,  slug-ftr 

K  a  constant,^  /2Tm  orya/2(Tm+  C) 

m  "jj     2f^/ITm,  subscript  meaning  maximum 

o  subscript  meaning  output 

R  resistance,  megohms  unless  otherwise  specified 

t  time,  seconds 

T  torque,  lb. -ft. 

v  voltage 

x  ^     I/2Tm    € 

£  error,  G.-90,  radians 

Cx  error  to  torque-saturate  the  servo,  Tm/u 

£  d€/dt 

£m  maximum  error  rate,  Tm/f 

€  d26/dt2 

j  damping  ratio 

6  angular  position,  radians 

/x  gain  constant 

1T  3. 14159 


time  constant,  I/f 

undamped  natural  frequency,  (a)^   -a/\ 
SX  input  angular  frequency 
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OPTIMUM  PREDICTION  FOR  A 
TORQUE  SATURATING  SERVO 

INTRODUCTION 

The  design  of  the  controller  for  a  servomechanism  is  based  to  a  large 
extent  on  the  particular  type  of  input  function  and  the  requirements  on  the  servo 
such  as  response  time,  accuracy,  and  amount  of  overshoot.    Linear  methods  of 
control  are  extensively  used  as  they  are  easy  to  implement  and  analyze.    Since 
linear  control  does  not  necessarily  represent  the  optimum  design  in  most  cases, 
considerable  work  has  been  done  recently  to  improve  performance  by  nonlinear 
techniques. 

One  such  technique  involves  use  of  prediction  or  switching  circuits  intended 
to  give  optimum  time  response  and  limitation  of  overshoot  for  step  or  ramp  in- 
puts or  combinations  of  the  two.    Previous  investigations  of  methods  of  control, 
using  some  measure  of  prediction  include  nonlinear  damping  functions  *    '    \ 
prediction  with  relay  servos*C'    '  and  multiple  mode  techniques*   \ 

Of  particular  interest  is  the  work  of  Hopkin'    ',  whose  approach  to  prediction 
is  very  similar  to  that  in  this  study  although  mechanized  in  a  different  manner. 

It  is  the  purpose  of  this  investigation  to  extend  previous  work  by  the 
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examination  of  use  of  a  particular  type  of  prediction  in  the  control  of  a  torque 
saturating,  second  order  servomechanism.    The  basic  theory  was  originally 
developed  by  L.  L.  Rauch  and  R.  M.  Howe*   '  and  is  summarized  in  Appendix  I. 
This  study  develops  a  method  of  implementation  of  the  theory  for  simulation  in 
several  phases  and  application  to  a  specific  instrument  servo. 

This  investigation  was  undertaken  by  William  W.  Gay  and  Wayne  S.  McCord 
as  a  joint  thesis  project  to  satisfy,  in  part,  the  requirements  for  a  Master  of 
Science  Degree  in  Aeronautical  Engineering  at  the  University  of  Michigan^ 

SIMULATION 

Simulation  was  undertaken  to  obtain  computer  solutions  to  theoretical 
equations  of  the  servo  problem  under  study.    These  solutions  were  obtained  for 
several  variations  of  the  servo  setup  in  order  to  make  a  series  of  comparisons. 

PROCEDURE 

The  servo  setups  considered  were  those  of: 

1 .  A  servo  using  proportional  and  error  rate  control  without  prediction. 

2.  A  servo  using  proportional  and  error  rate  control  with  prediction. 

3.  The  same  servo  as  in  #2  but  with  an  approximate  prediction  curve. 

a.  No  friction. 

b.  With  Coulomb  and/or  viscous  friction 

(but  no  compensation  for  friction  in  the  prediction  curve). 

c.  With  Coulomb  and/or  viscous  friction 

(with  compensation  for  Coulomb  or  viscous  friction). 

4.  The  same  servo  as  in  #3,  a. ,  but  with  approximate  differentiation  of 
error  to  obtain  error  rate. 


The  first  point  considered  in  the  simulation  was  to  show  the  improvement 
in  response  time  and  reduction  of  overshoot  by  adding  prediction  to  a  servo  with 
proportional  plus  error  rate  control. 

The  equation  for  a  torque-saturating,  frictionless  servo  employing  proportional 
and  error  rate  control  is 


le.^.QC^ 


To  adapt  the  equation  to  the  computer  for  solution,  it  was  time  scaled  by  the 
undamped  natural  frequency,  a)   ,  and    magnitude  scaled  by  the  error  required 

for  torque  saturation,  €.#  .    These  quantities  are  defined  as 

60  2      =      J±- 
n  s\ 

T 
€,    :    J^ 

*       A 

Also,   Ce  r  U. 

where  y  is  the  damping  ratio.    The  equation  can  then  be  written 


The  arbitrary  values  of  parameters  chosen  were  con  =   1,    ^=0.5,  and 
€.x  =   0. 025,  where  unity  was  represented  by  100  volts  in  the  computer.    The 
solutions  were  then  presented  in  dimensionless  form  as  €/con  versus  eMA  in 
the  phase  plane  and  e/e^  versus  oont   in  the  time  plots.    Figure  1  shows  the 
computer  circuit  for  solution  of  this  equation. 

The  effect  of  torque  saturation  was  obtained  by  limiting  with  operational 
amplifiers.    The  output  stage  of  each  amplifier  was  a  cathode  follower.    When 


4 
the  outputs  of  two  amplifiers  were  connected  in  parallel,  the  more  positive  signal 
determined  the  output  of  the  pair.    Thus  by  introducing  a  constant  positive  voltage 
into  the  input  of  one  of  the  paralleled  amplifiers,  the  output  of  the  pair  was  limited 
in  the  negative  direction.    By  using  two  such  pairs  in  cascade,  both  plus  and  minus 
limiting  was  obtained. 

The  addition  of  prediction  modifies  the  equation  of  motion  of  the  above  servo  to 
the  form  developed  in  Appendix  I. 


re.^e+Cc^at^+^N^^i 


The  method  for  adapting  this  equation  to  a  suitable  form  for  computer  solution  is 
identical  to  the  procedure  used  previously  with  the  result 
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This  equation  was  actually  solved  by  two  different  computer  circuits.    Figure 
2  shows  the  computer  circuit  for  exact  solution.    A  servo  multiplier  was  used  to 
form  the  product  I—  I  £-. 

The  computer  circuit  in  Figure  3  was  also  used  to  solve  this  equation.    In 
this  case,  a  diode  limiter  network  at  the  input  to  an  operational  amplifier  generated 
an  approximation  to  the  prediction  curve.    The  limiting  network  and  approximate 
prediction  curve  for  this  case  with  no  friction  are  shown  in  Figure  5.    The  diode 
networks  used  in  this  and  subsequent  computer  circuits  were  developed  by  the 
methods  presented  in  Appendix  n. 

Friction  was  introduced  into  the  servo  problem  with  no  change  in  the  pre- 
diction curve. 

With  Coulomb  or  viscous  friction,  the  equations,  when  adapted  for  computer 


solution,  became  respectively 


4  4^4>^+£MU 

where  *r  =   I/f  is  the  servo  time  constant,  and  C  is  the  Coulomb  friction  torque. 

The  effect  of  Coulomb  friction  was  introduced  in  the  computer  as  a  change 
in  torque  saturation  levels  depending  on  the  sign  of    9    .    Viscous  friction  was 
introduced  in  the  computer  as  a  feedback  in  the    6 ^ 'oo     integrator.    The  circuit 
additions  are  shown  in  Figure  3  in  dashed  lines. 

Next,  compensation  for  Coulomb  and  viscous  friction  effects  was  introduced 
by  changing  the  approximate  prediction  function. 

When  the  prediction  is  optimized  for  a  specific  Coulomb  friction,  the  equation 
for  the  predictor  servo,  when  adapted  for  computer  solution  becomes 


ft  -  [« *  «"&*,«,+  ztSratel^L,,  * 


The  diode  limiting  network  and  approximate  prediction  curve  designed  to  include 
compensation  for  Coulomb  friction  of  C/T     =    0.  3  were  unchanged  from  Figure  5. 
The  only  requirement  was  to  multiply  the  coefficient  of  k/u>    by  the  factor  t/  *  •  Z 

When  the  prediction  curve  is  optimized  for  a  specific  viscous  friction,  the 
equation,  when  adapted  for  computer  solution  becomes 


The  diode  limiting  networks  and  approximate  prediction  curve  compensating  for 
viscous  friction  are  shown  in  Figure  6  for  T  u)    =   10  and  T  o>  n  =   8.  95. 

Approximate  differentiation  was  introduced  only  for  the  servo  with  no  friction. 
The  computer  equation  developed  previously  still  applies.    The  computer  setup 
and  approximate  differentiation  circuit  are  shown  in  Figure  4. 

RESULTS  AND  DISCUSSION 

The  results  of  the  simulation  study  were  derived  entirely  from  the  phase  plane 
and  time  plots  for  the  various  conditions  outlined  previously.    Since  some  effects 
are  more  apparent  in  one  representation  than  the  other,  certain  solutions  were 
recorded  in  the  phase  plane  and  also  as  a  time  function.    Complete  duplication  of 
phase  plane  and  time  plots  was  not  attempted-    The  solutions  were  grouped  so  a 
comparison  of  certain  methods  of  control  could  be  made  to  determine  the  answers 
to  five  questions. 

1.  Does  the  use  of  this  type  of  prediction  materially  improve  the  performance 
of  a  servo  system? 

2.  Does  an  approximation  to  the  prediction  curve  materially  degrade  servo 
performance? 

3.  How  much  does  friction  affect  servo  performance  when  prediction  is  used? 
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4.  Does  compensation  for  friction  effects  materially  improve  the  servo 
performance? 

5.  Is  approximate  differentiation  of  error  to  obtain  error  rate  feasible  for 
this  system? 

A  measure  of  the  improvement  in  performance  due  to  prediction  is  available 
from  the  solutions  in  Figures  7  through  13.    For  small  step  inputs  the  servo  with 
prediction  shows  no  reduction  in  the  size  of  the  overshoot  or  response  time. 
However,  as  the  size  of  the  input  was  increased,  the  improvement  became  very 
pronounced.    For  a  step  of  e  /e^     =   32,  the  overshoot  was  reduced  by  about 
95%  and  response  time  by  65%. 

Of  interest  in  this  comparison  are  the  responses,  in  Figures  14  through  16, 
of  the  servo  with  and  without  prediction  to  sine  wave  inputs.    Although  there  is  no 
particular  significance  to  this  type  of  response  and  the  prediction  was  not  optimized 
for  it,  the  improvement  was  still  evident.    For  the  highest  frequency  input,  the 
servo  with  prediction  could  follow  the  signal,  while  the  servo  without  prediction 
developed  a  divergent  error.    This  condition  was  shown  in  Figure  16  where  the 
input  frequency  was  1/3  of  the  servo  natural  frequency  and  the  ratio  of  the  maximum 
acceleration  of  the  input  to  the  maximum  acceleration  of  the  servo  was  8/9. 

Typical  solutions  of  the  servo  response  using  an  approximation  to  the  pre- 
diction curve  are  presented  in  Figures  17  (a)  and  (b)  and  18.    A  comparison  of  these 
solutions  to  those  in  Figures  9  (a),  8  (a)  and  12  respectively  shows  no  preceptible 
difference.    Since  the  operation  of  the  system  with  this  five-segment  approximation 
was  satisfactory,  no  solutions  were  made  using  a  more  exact  approximation  to  the 
true  switching  curve.    All  the  remaining  solutions  were  obtained  with  this  method 
of  generating  the  prediction  curve. 


8 

The  selection  of  friction  coefficients  to  demonstrate  the  effects  of  friction 
was  somewhat  arbitrary.    The  Coulomb  friction  to  torque  ratio  corresponds  to 
that  specified  for  the  servo  later  tested.    Since  the  viscous  friction  of  this 
system  was  small,  corresponding  to  T  to     -   100,  additional  values  of  viscous 
friction  were  selected  to  emphasize  the  effect.    The  change  in  the  response  with 
either  type  of  friction  is  most  evident  in  the  phase  plane  plots  of  Figures  19  and 
20.    This  change  appears  as  a  considerable  reduction  in  error  rate,  but  it  may 
be  observed  in  the  error  versus  time  plots  in  Figures  21  and  22  that  the  actual 
increase  in  response  time  was  small  even  for  these  large  amounts  of  friction. 
It  may  be  noted  that  when  either  type  of  friction  was  present  the  solution  path 
lay  somewhat  within  the  linear  region  during  deceleration,  indicating  that  less 
than  maximum  torque  was  being  used. 

In  view  of  the  fact  that  friction  has  the  effect  of  slowing  servo  response  to 
some  extent,  it  was  necessary  to  determine  if  response  could  be  improved  by 
changing  the  prediction  curve  to  permit  utilization  of  maximum  torque  during 
deceleration. 

Figures  23  and  24  show  phase  plane  and  time  plots  where  the  prediction 
curve  was  altered  to  take  into  account  the  presence  of  Coulomb  friction.    These 
plots  may  be  compared  to  those  in  Figures  20  and  22  for  the  same  conditions 
but  with  prediction  calculated  for  no  friction.    A  very  minor  improvement  in 
response  time  was  obtained  by  altering  the  prediction  curve  to  compensate  for 
Coulomb  friction.    More  improvement  was  noted  as  size  of  the  step  was  increased. 

In  considering  the  various  runs  in  Figures  25  through  29  it  immediately  be- 
comes obvious  that  little  improvement  in  servo  time  response  was  obtained  when 
the  prediction  curve  was  altered  to  take  into  account  viscous  friction.    It  may  be 
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noted  that  Figures  25  and  26  show  a  case  where  no  velocity  saturation  due  to 
viscous  friction  occured.    Figures  27  through  29  show  a  case  where  velocity 
saturation  was  approached. 

It  may  be  said  that  prediction  to  include  either  type  of  friction  effects  is 
feasible  but  offers  little  advantage. 

By  examining  Figure  30  it  is  apparent  that  approximate  differentiation  was 
satisfactory  in  that  the  phase  plane  plot  appears  quite  similar  to  phase  plane 
plots  where  the  true  value  for  error  rate  was  used  in  computation. 

Stability  was  a  problem  in  the  simulation,  however,  and  would  presumably 
be  a  serious  problem  in  application.    The  time  constants  of  the  differentiating 
circuit  for  the  simulation  were  selected  in  order  to  avoid  oscillation. 

APPLICATION  TO  AN  INSTRUMENT  SERVOMECHANISM 

A  physical  system  was  constructed  to  demonstrate  in  practice  the  use  of 
approximate  differentiation  and  an  approximate  prediction  curve  similar  to  the 
one  examined  in  the  simulation. 


EQUIPMENT 

Motor.    The  servo  motor  was  a  two  phase,  400  cycle  per  second,  115  volt 

Transicoil  motor,  type  18M.    Rated  stall  torque  -3.0  oz.  -in. ,  moment  of  inertia 

2 
0. 0197  oz.  -in.  ,  stall  watts  per  phase  -  18  watts,  Coulomb  friction  -  0. 10  oz.  -in. 

coefficient  of  viscous  friction  -  0.  00285  oz.  -in.  -sec. 

Potentiometer.  Feedback  voltage  was  obtained  from  a  three  section,  con- 
tinuous film  potentiometer,  model  number  205,  made  by  Computer  Instruments 
Corporation,  Resistance  -  50, 000  ohms  per  section.    Total  moment  of  inertia  - 
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2 
0. 164  oz.  -in.    ,  total  Coulomb  friction  -  0.  895  oz.  -in. ,  total  viscous  friction  - 

negligible. 

Diode  limiting  circuits.    This  is  the  circuit  described  in  Appendix  n, 
Figure  n-3. 

Compensation  and  control  circuits.    The  compensation  and  control  circuits 
utilized  the  operational  amplifiers  of  a  standard  analog  computer. 

Power  Amplifier.  Two  different  power  amplifiers  were  used  to  drive  the 
control  phase  of  the  motor.  One  was  a  115  volt,  400  cycle  per  second,  push- 
pull,  magnetic  amplifier,  Magnetic  Amplifiers  Inc.    Model  MA  41501-CY. 

The  other  amplifier  was  an  A.  C.  electronic  amplifier  used  in  conjunction 
with  a  400  cycle  per  second  chopper  to  convert  the  D.  C.  control  signal  to  a 
400  cycle  per  second  A.  C.  signal.    This  amplifier  was  originally  designed  as  a 
60  cycle  per  second  servo  amplifier  and  was  adapted  for  use  with  this  system. 

The  choice  of  equipment  was  governed  primarily  by  the  requirements  of  the 
control  circuits,  specifically  the  differentiator.    The  film  potentiometer  was 
selected  in  preference  to  a  wire  wound  potentiometer  because  the  signal  obtained 
had  less  noise  and  better  resolution.    The  wire-to-wire  discontinuities  of  a  wire 
wound  potentiometer  would  have  introduced  large  errors  in  the  differentiated 
signal  and  contributed  to  the  instability  of  the  differentiating  amplifier. 

The  motor  was  selected  from  those  available  on  the  basis  of  largest  torque 
to  inertia  ratio.    It  was  coupled  directly  to  the  potentiometer  to  eliminate  the 
backlash  and  binding  problems  of  gear  trains.    Although  the  motor-to-load 
inertia  balance  was  not  good,  the  acceleration  loss  due  to  direct  coupling  was 
small.    The  best  gear  ratio,  if  used,  would  have  been  about  three  to  one. 

Initially  a  magnetic  amplifier  was  chosen  for  power  amplification  because 
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of  simplicity  and  reliability.    The  severe  bandwidth  limitations  were  alleviated 
somewhat  by  compensation  as  described  in  Appendix  III  so  that  the  response  was 
flat  out  to  200  cycles  per  second.    However,  the  response  time  of  the  servo  was 
short  enough  that  the  200  cycles  per  second  bandwidth  still  represented  a  limi- 
tation which  appeared  as  a  time  delay.    When  the  system  gain  was  increased  to 
obtain  reasonable  resolution,  the  time  lag  in  the  magnetic  amplifier  sustained 
oscillations  of  the  servo  at  a  frequency  of  80  to  120  cycles  per  second,  depending 
upon  the  particular  gain  setting.    The  oscillations  were  similar  to  those  that 
would  be  experienced  by  an  identical  bang -bang  servo  with  a  fixed  switching 
delay  of  about  5  milliseconds.    This  time  delay  was  verified  by  observation  of 
the  control  and  output  signals  of  the  magnetic  amplifier  on  a  dual  beam  oscillo- 
scope.   The  electronic  amplifier  and  chopper  were  substituted  when  operation 
with  the  magnetic  amplifier  did  not  produce  the  potential  resolution  or  tracking 
performance  of  the  servo. 

PROCEDURE 

The  servo  was  operated  to  obtain  a  comparison  of  the  system  with  pro- 
portional control,  proportional  plus  error  rate  control  and  proportional  plus 
error  rate  control  with  prediction.    The  comparison  was  made  by  observing 
only  the  step  response,  and  no  attempt  was  made  to  analyze  the  servo  completely. 
As  an  item  of  interest,  the  velocity  lag  and  constant  input  errors  were  measured. 

The  basic  servo  circuit  diagram  is  shown  in  Figure  31.    Appropriate  modi- 
fications to  this  basic  circuit  were  made  to  obtain  the  particular  types  of  control. 
The  servo  loop  gain  was  adjusted  to  the  highest  possible  value  consistent  with 
stable  operation.    If  the  gain  was  further  increased,  the  servo  went  into  a  very 
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low  amplitude  oscillation  at  about  160  cycles  per  second. 

The  time  constants  of  the  approximate  differentiating  circuit  were  determined 
empirically  with  a  view  toward  limiting  noise  and  high  frequency  oscillations  yet 
retaining  sufficient  bandwidth  for  good  operation.    The  time  constant  of  0.  0005 
seconds  corresponds  to  319  cycles  or  roughly  twice  the  maximum  response 
frequency  of  the  servo.    The  basic  differentiating  circuit,  with  constant  gain  at 
high  frequencies  (within  the  capabilities  of  the  operational  amplifier),  passed 
too  much  noise,  so  the  additional  capacitor  was  added  in  the  feedback  to  give  a 
-6  db  per  octave  roll  off  of  gain  for  frequencies  above  319  cycles  per  second. 

The  approximations  in  differentiation  modified  the  error  rate  signal  such 
that  the  amplitude  was  not  increased  in  proportion  to  the  frequency  as  prescribed 
for  exact  differentiation.    Some  compensation  for  the  resulting  loss  of  gain  was 

made  in  the  next  amplification  stage  where  the  product  _  /  _1T    *  e    was 

V2<Tm-^C) 

formed.    In  place  of  using  a  calculated  gain  of  3. 16  in  this  amplifier,  the  gain 
was  increased  to  4.    The  value  of  4  was  determined  empirically  for  satisfactory 

servo  response.    The  factor  rf   in  -^      ,TT„  1  is  the  scale  factor  between  the 

V     2(Tm+  C) 

output  measured  in  radians  and  computer  units,  where  100  volts  represents  unity 
in  the  computer. 

Photographs  of  the  cathode  ray  oscilloscope  presentation  were  made  for 
phase  plane  and  error  versus  time  response  of  the  servo  to  a  100  volt  peak-to- 
peak  3  cycle  per  second  square  wave  input.    Similar  photographs  were  made  of 
the  error  versus  time  response  for  a  10  volt  peak-to-peak  5  cycle  per  second 
square  wave  input. 

The  constant  input  error  was  measured  by  disturbing  the  servo  slightly 
and  noting  the  maximum  error  when  the  output  had  come  to  rest.    The  error 
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was  actually  measured  on  a  voltmeter  having  a  full  scale  deflection  of  0. 1  volt. 

The  velocity  lag  was  measured  in  a  manner  similar  to  the  constant  input 
error.    The  effect  of  a  ramp  input  was  produced  by  a  0.  5  cycle  per  second 
triangular  wave.    The  frequency  of  the  triangular  wave  was  low  enough  that  the 
transient  effects  damped  sufficiently  to  permit  accurate  measurement  of  the 
velocity  lag. 

RESULTS  AND  DISCUSSION 

Qualitatively  the  performance  of  the  servo  corresponded  closely  to  results 
of  the  simulation.    Although  specific  performance  of  the  servo  was  not  an 
objective,  the  response  time,  ability  to  track  a  slowly  varying  signal  and  the 
resolution  approached  closely  the  potential  capability  of  the  system. 

The  measured  value  of  maximum  steady  state  error  was  ±  0.  005  volts. 
The  terminal  voltages  on  the  potentiometer  were  ±  100  volts  with  the  center- 
tap  grounded.    Based  on  this  reference  voltage  range  of  200  volts,  the  maximum 
steady  state  error  corresponds  to  resolution  of  at  least  one  part  in  40,  000 

From  the  maximum  steady  state  error,  an  approximation  to  the  width  of 
the  linear  region  (2  £^  )  may  be  calculated.    The  Coulomb  friction  to  rated 
torque  ratio  was  0.  33.    However,  actual  measurement  of  the  stall  torque  of  the 
motor  as  driven  by  the  electronic  amplifier  was  about  2  oz.  -in. ,  giving  an  actual 
friction  to  torque  ratio  of  about  0. 5.    Using  this  ratio,  the  linear  region  was 
calculated  to  be  in  the  order  of  0.  02  volts. 

Among  the  factors  contributing  to  the  low  value  of  torque  were  poor  output 
waveform  from  the  electronic  amplifier  and  somewhat  less  than  90°  phase  shift 
between  the  two  motor  phase  voltages. 
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The  photographs  of  the  oscilloscope  presentations  of  the  phase  plane  and 
error  versus  time  responses  are  shown  in  Figure  32.    The  system  using  propor- 
tional control  exhibited  a  series  of  fairly  large  overshoots  and  long  response 
time.    Addition  of  error  rate  control  (  J  =  0.  5)  reduced  the  size  and  number  of 
overshoots  and  response  time.    Furthermore,  addition  of  prediction  essentially 
eliminated  overshoot  and  very  substantially  reduced  response  time.    As  expected, 
the  improvement  due  to  prediction  was  much  more  pronounced  for  the  larger 
steps.    As  an  example,  for  the  100  volt  step  input  the  response  times  were  0. 17, 
0. 13  and  0.  06  seconds  for  the  three  types  of  control. 

Measurements  were  made  of  velocity  lag  and  are  recorded  in  Table  I. 


TABLE  I 


VELOCITY  LAG  (All  Types  of  Control) 


0    (volts/sec. ) 

Error   (volts) 

0.5 

0.01 

1 

0.014 

2 

0.02 

5 

0.03 

10 

0.04 

20 

0.06 

40 

0.076 

80 

0.08 

120 

0.10 

160 

0.13 

Since  the  rotation  of  the  motor  was  limited  to  at  most  one  revolution,  the 
velocity  never  exceeded  20%  of  the  maximum  motor  velocity.    Within  this  range, 
the  viscous  friction  effects  of  the  servo  were  negligible.    For  practical  con- 
siderations, this  fact  plus  the  very  narrow  linear  region  makes  the  response  of 
the  servo  to  a  step  input  of  several  volts  very  similar  to  that  of  a  bang-bang  servo 
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with  only  Coulomb  damping. 

The  response  of  the  servo  with  proportional  control  to  the  100  volt  step 
input  exhibits  a  ratio  of  successive  overshoots  of  about  0.  4.    Calculations 
using  this  ratio  for  a  theoretical  bang-bang  servo  with  Coulomb  damping 
indicate  a  Coulomb  friction  to  torque  ratio  of  0. 43,  which  checks  closely  with 
the  ratio  of  the  measured  values,  0.5. 

The  response  of  the  servo  with  prediction  can  also  be  compared  to  the 
theoretical  response  of  a  corresponding  bang-bang  servo.    Again  for  the  100 
volt  step  input,  the  observed  response  time  was  0. 06  seconds  which  compares 
to  a  computed  time  of  0.  063  seconds. 

CONCLUSIONS 

Based  on  the  results  obtained  in  the  simulation  and  in  application  to  the 
instrument  servo  it  is  concluded  that: 

1.  Use  of  prediction  reduces  response  time  and  overshoot,  particularly 
for  large  step  and  ramp  inputs. 

2.  Use  of  prediction  improves  transient  response  of  the  servo  for  arbitrary 
input  signals,  such  as  a  sine  wave  input,  even  though  the  prediction  was 
optimized  only  for  step  and  ramp  inputs. 

3.  Use  of  simple  approximations  to  the  exact  prediction  curve  does  not 
degrade  performance. 

4.  Modification  of  the  prediction  curve  to  optimize  for  friction  effects 
was  possible  but  little  improvement  was  achieved  over  the  prediction 
curve  optimized  for  no  friction 

5.  Noise  level  in  the  output  signal  from  the  film  potentiometer  was  low, 
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thus  permitting  differentiation  with  an  operational  amplifier  to  obtain 
approximate  error  rate  signal. 

6.  Resolution  of  the  film  potentiometer  and  direct  coupling  permitted  a 
degree  of  control  sufficient  to  give  servo  resolution  of  one  part  in 
40, 000. 

7.  For  high  loop  gains  the  region  of  linear  operation  was  very  narrow, 
hence,  the  system  acted  much  like  a  bang-bang  servo  with  Coulomb 
damping. 
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FIGURE  5  —  (a)    Approximate  Prediction  Curve  for  a  Frictionless 
Servo,  K  =  yu/2Tm,  and  a  Servo  with  Coulomb 
Friction,  K  =yx/2(Tm  f  C) 
(b)     Corresponding  Diode  Function  Generator  Circuit 
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FIGURE  6  —  (a)  Approximate  Prediction  Curve  for  a  Servo  with 
Viscous  Friction,  tu3q=  10,    8.95 
(b)    Corresponding  Diode  Function  Generator  Circuit 
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^IGURE  7  —  Step  Response  of  a  Frictionless  Servo  with 

Proportional  and  Error  Rate  Control,  C     -  lA<->n 


(a)  With  Prediction 


(b)  Without  Prediction 
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FIGURE  8  —  Step  Response  of  a  Frictionless  Servo  with  Proportional 
and  Error  Rate  Control,  Ce  -    l/o> 
(a)   With  Prediction  (b)   Without  Prediction 
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FIGURE   9  —  Step  Response  of  a  Frictionless  Servo  with  Proportional 
and  Error  Rate  Control,  Ca    -   l/oon 


(a)  With  Prediction 


(b)  Without  Prediction 
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FIGURE  10  —  Step  Response  of  a  Frictionless  Servo  with 

Proportional  and  Error  Rate  Control,  Ce    -   l/6n 


(a)    With  Prediction 


(b)    Without  Prediction 
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FIGURE  11  —  Ramp  Response  of  a  Frictionless  Servo  with 

Proportional  and  Error  Rate  Control,  Ce  =  l/u>  n, 


e./ojQ  s  o.i 

(a)   With  Prediction 


(b)   Without  Prediction 
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FIGURE  14  —  Sine  Wave  Response  of  a  Frictionless  Servo  with 
Proportional  and  Error  Rate  Control,  Ce  =   lA>n, 
Jv^n:l/5,    e'i  (max)/<9o(max)    =    0.32 
(a)    With  Prediction  (b)  Without  Prediction 
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FIGURE  15  —  Sine  Wave  Response  of  a  Frictionless  Servo  with 
Proportional  and  Error  Rate  Control,  Ce  =   \/cO  n, 
-0/wn  -  1/4,    6i  (max)/  6>0(max)  =  1/2 
(a)    With  Prediction  (b)    Without  Prediction 
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FIGURE  16  —  Sine  Wave  Response  of  a  Frictionless  Servo  with 
Proportional  and  Error  Rate  Control,  Ce  =    1/a)  n, 
-%^n  =  1/3,    '6i  (max)  /  9Q  (max)  =  8/9 
(a)  With  Prediction  (b)    Without  Prediction 
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FIGURE  17  —  Step  Responses  of  a  Frictionless  Servo  with 

Proportional  and  Error  Rate  Control,  Ce  =  l/60n, 
Plus  Approximate  Prediction. 
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FIGURE  21  —  Step  Response  of  a  Servo  with  Prediction  Optimized 
for  No  Friction 
(a)   No  Friction      (b)  fto     -  100     (c)  foJ     =  10 
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FIGURE  22  —  Step  Response  of  a  Servo  with  Prediction  Optimized 
for  No  Friction     (a)    No  Friction      (b)   No  Viscous 


Friction,    C/T 


0.3     (c)ro>n  =  l0,    C/Tm  =  0.3 
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FIGURE  24  —  Step  Response  of  a  Servo  with  Prediction 

Optimized  for  Coulomb  Friction,  C/T  m  =  0. 3 


(a)    No  Viscous  Friction 
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FIGURE  26  —  Step  Responses  of  a  Servo,  Too       =   10,  with  Prediction 
Optimized  for  n 

(a)   No  Friction         (b)  Viscous  Friction,  Too      -    10 
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FIGURE  29  —  Step  Responses  of  a  Servo  Approaching  Velocity 

Saturation,  Ta)     =8.95,    Prediction  Optimized  for 

(a)   No  Friction     (b)    Viscous  Friction,  red       -   8.95 
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APPENDIX  I 
OPTIMUM  SWITCHING  CRITERIA 

The  minimum  response  time  for  a  servo  system  with  limited  torque  is 
obtained  if  the  maximum  available  torque  is  applied  to  accelerate  the  system, 
and  then  the  torque  reversed  for  maximum  deceleration  to  arrive  simultaneously 
at  zero  error  and  zero  error  rate.    If  the  torque  is  then  removed,  the  error  will 
remain  zero  until  another  command  signal  is  given  to  the  system.    This  type  of 
response  is  discussed  for  three  categories;  a  system  with  no  friction,  a  system 
with  only  Coulomb  friction,  and  a  system  with  only  viscous  friction. 

The  problem  of  determining  the  exact  time  to  reverse  the  torque  depends 
in  part  upon  the  characteristics  of  the  command  signal.    For  the  first  two  cate- 
gories the  acceleration  of  the  command  signal  is  assumed  to  be  zero.    That  is, 
the  switching  time  is  optimized  for  constant  or  uniformly  varying  command 
signals.    For  the  case  with  viscous  friction,  the  switching  time  is  optimized 
for  a  constant  command  signal. 

First  consider  the  system  with  no  friction.    The  equation  of  motion  is 

I^o  =  ±Tm=   tonJue 
or,  in  terms  of  error,  since   0^  =   0, 

i 

Integrating  two  times  and  assuming  initial  conditions  of  e  (0)   =    0  and 
€  (0)   =   0,  the  solution  of  the  equation  is 

2 

c        -     Tm         t 

£z  +  "i T 
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or 


2T 


since  t 


€     from  the  first  integration. 


In  the  phase  plane  of  error  and  error  rate,  the  solution  paths  are  two 
parabolas  passing  through  the  origin  as  shown  in  Figure  1-1. 


FIGURE  1-1 
Phase  Plane  Plot  of  Solutions  Through  the  Origin 


The  parts  of  interest  lie  in  the  second  and  fourth  quadrants,  the  other  two  portions 
represent  divergent  solutions.    The  phase  plane  is  thus  divided  into  two  regions 
by  the  two  segments  of  parabolas  comprising  the  switching  curve.    The  equation 
for  this  curve  is 


£  = 


I 
2TT 


■|€|€ 
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For  a  solution  starting  above  and  to  the  right  of  the  switching  curve,  full  positive 
torque  applied  will  cause  the  solution  to  move  along  a  parabolic  path  until  it 
reaches  the  switching  curve  in  the  fourth  quadrant.    If  the  torque  is  then  reversed, 
the  response  will  follow  the  switching  curve  to  the  origin,  and  the  torque  can  be 
reduced  to  zero.    The  same  procedure  applies  to  the  other  half  of  the  phase  plane 
with  the  opposite  sense  of  torque.    Sample  phase  plane  trajectories  are  shown  in 
Figure  1-2. 


FIGURE  1-2 
Example  Trajectories 


Although  this  switching  criterion  represents  the  fastest  response  for  large 
errors  or  error  rates,  a  practical  system  mechanized  in  this  manner  would  not 
follow  slowly  varying  inputs  smoothly  due  to  switching  delays  and  other  non-ideal 
characteristics.    A  linear  system  with  error  rate  or  derivative  feedback  in 
addition  to  proportional  control  is  better  for  this  type  of  operation  near  the  origin 
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in  the  phase  plane.    The  torque  in  a  linear  system  with  proportional  plus  error 
rate  control  is  proportional  to    £  +-C    e   ,  and  the  optimum  nonlinear  system 
utilizes  maximum  torque  tending  to  satisfy  the  equation 

£  +  — —  |6|€       -  0 
2Tm 

In  combining  the  two  systems,  the  torque  should  be  made  proportional  to 

C  +  C   €  +  _L_|£|€  .    Near  the  origin  of  the  phase  plane,  €  is  small 
©  2T 

and  the  term_L_  je|  £.       is  negligible  compared  to  the  others,  hence  the  servo 

2Tm 
response  is  linear  as  desired.    Away  from  the  origin  the  CQ€   term  must  be 

limited  so  the  response  will  tend  to  make 

For  a  large  gain  constant,  /u  ,  a  small  error,    €     =   TmAoc  ,  will  torque 
saturate  the  system.    If  the  Ce<£     term  is  made  to  saturate  at  a  level  exactly 
equal  to  the  error  required  to  torque  saturate  the  servo,  then  away  from  the  origin 
the  response  will  lie  exactly  on  the  optimum  switching  curve.    The  equation  of 
motion  of  the  servo  is  then 


10.  V^+iC.eU,^  ^lelel      T 


where  sat  ±  T    /m    represents  the  limits  on  the  particular  quantity  in  brackets. 
In  the  phase  plane,  the  torque  equation  represents  a  switching  curve  with  a 
linear  region  adjacent  to  the  switching  curve.    The  limits  of  this  linear  region 
are  shown  in  Figure  1-3.    The  zero  torque  line  lies  midway  between  the  boundaries 
of  the  linear  region  in  the  no  friction  case. 
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FIGURE  1-3 
Phase  Plane  Trajectory  For  the  Optimum 
Switching  Curve  with  Linear  Region 


The  extension  of  the  above  development  to  the  case  of  a  system  with  Coulomb 
friction  is  quite  straightforward.    The  Coulomb  friction  represents  a  constant 
torque,  C,  which  adds  to  servo  torque  during  deceleration  and  subtracts  from 
servo  torque  during  acceleration.    The  switching  curve  equation  then  becomes 


e  +• 


■le|€       =.o 


2<Tm-hC)' 
and  the  corresponding  equation  of  motion  of  the  system  is 


/" 


r 


The  effects  of  Coulomb  friction  on  the  switching  curve  are  to  increase  the 
magnitude  of  the  slope  in  a  ratio  (1  +-  C/Tm),  and  to  shift  the  zero  torque  line 
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from  the  center  of  the  linear  region. 

In  developing  the  equation  of  the  switching  curve  for  a  system  with  viscous 
friction,  the  additional  condition  of  6  .    -   0  was  applied.    The  equation  of  motion 
of  the  system  is  now 


ieQ+    f9o   =   torque   =  ±  Tm 


or  in  terms  of  error 


r  £  -f-  £  =  r  e 

+      m 


where  r  =  I/f  =   time  constant  and  £      =   T    /f  =   maximum  error  rate 


=   T^ml1-^7]    ■     e<0)   =    0 


or      t  =   -rln(l  t  £/em) 


£  =    +  ^  m  [*  "  r  <X  "  e"^r    >]  •  €  (0)   =    0 

6=  ±  dmr[in<i*€/€m>-  e/em] 


e+ri  =remin(i^  |e/ej>.  €>o 


e+re  =  -t£     i*<i  +  |e/e J),  £<o 

m  m1 


The  equation  of  motion  of  the  servo  system  is  then 
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To  compare  the  switching  curve  for  viscous  friction  to  the  one  with  no 
friction,  consider  for  the  moment,  positive  G, . 
Let*  =    iV2Tm       £ 
andm  =   ^2f2/lTm   =     p/T  €  m 
The  switching  curve  for  no  friction  is 

I        -2  2 

fc       =   -  x 


2T 
*xm 


The  switching  curve  for  viscous  friction  is 

£-   2/m2  [ln(l  +  mx)   -   mx] 
Since  mx  <  1,  expand  ln(l  -t-  mx)  in  a  series 

e=-i«{[— -^2+<r3-<^4+ ]  --} 

2        r     W^  (m«)3       <W  ,  1 

6:  ^2  L'.T  +-s r-  + J 

e=-^    [1.^L+i"l'.   W+.....] 


It  is  apparent  that  for  no  friction,  m  «   0,  the  two  curves  are  identical. 
The  difference  is  important  only  when  mx  is  significant  compared  to  1.    The  same 
result  can  be  obtained  by  comparing  the  slopes  of  the  two  curves. 

4^-=   -2x  forf  =   0 

and*" 

de  -2x 

&rz     1  -h  mx        ^  f  ^  0 
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Comparative  curves  are  shown  in  Figure  1-4  for  several  values  of  m. 
The  curves  are  terminated  at  mx  ^  1. 
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APPENDIX  H 
THE  FUNCTION  GENERATOR 

Any  single  valued  function  can  be  approximated  by  straight  line  segments 
to  a  reasonable  degree  of  accuracy  with  an  operational  amplifier  and  biased 
diodes.    The  problem  is  only  one  of  calculating  the  resistance  values  and  bias 
voltages  required  to  give  the  proper  slope  and  location  of  the  straight  line 
segments. 

Consider  the  circuit  and  response  curve  in  Figure  II- 1,    Make  the  assumption 
that  the  resistances  are  large  compared  to  the  grid  to  plate  resistance  of  the 
diode  and  the  power  supply  impedance  so  that  the  voltage  drop  across  the  con- 
ducting diode  and  in  the  power  supply  are  negligible. 


eo 


*1 
-R 


(-v) 


ei>  ' 


Rl    +R2 


Ri   +  R2 


Rt     +   R2  1 

FIGURE  n-1 
(a)   Diode  Limiting  Circuit  (b)   Response  Curve 


If  another  input  is  now  added  with  the  diode  reversed  and  biased  positively, 
the  output  can  be  made  symmetrical  with  respect  to  the  origin.    Symmetry  is 
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not  necessary,  but  is  demonstrated  below  in  Figure  n-2  because  of  the  par- 
ticular function  of  interest  in  this  report. 


A  eL        Rz  K, 


=     0 


Rl   +   R2 

eii    <  R,  v 


eil    > 


Rl  +h2 


FIGURE  H-2 


(a)  Symmetrical  Diode 
Limiting  Circuit 


(b)  Response  Curve 


By  adding  other  inputs  with  different  resistances  and  bias  voltages,  more 
straight  line  segments  can  be  generated  and  combined  with  the  existing  ones. 
Each  input  contributes  a  change  in  slope  in  the  region  where  the  diode  is  cut 
off,  and  a  constant  output  elsewhere.    Thus,  by  choosing  the  proper  cut  off 
voltage  and  resistance,  the  point  of  slope  change  and  the  amount  of  change 
are  determined. 

For  most  problems  it  is  desirable  to  have  a  well  defined,  steady  zero 
reference.  If  the  diodes  are  conducting  for  zero  input,  drifting  of  the  bias 
voltage,  changes  in  the  diode  resistance  and  inequalities  of  resistances  tend 
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to  give  output  drift  and  balancing  problems.    For  these  reasons  it  is  preferable 
to  operate  with  the  diodes  cut  off  at  zero  input.    This  condition  in  general  re- 
quires an  additional  operational  amplifier  if  the  second  derivative  of  the  function 
is  positive,  as  it  is  in  this  servo  problem.    The  procedure  will  be  demonstrated 
for  generation  of  the  torque  switching  curve  used  for  the  frictionless  servo. 

It  should  be  noted  at  this  point  that  biased  diodes  could  be  used  in  the 
amplifier  feedback  circuit  with  an  effect  opposite  to  those  in  the  input.    This 
procedure  was  not  used  because  it  is  difficult  to  generate  a  zero  slope,  and  also 
a  change  of  slope  of  any  segment  involves  solution  of  a  set  of  "n"  simultaneous 
equations  to  determine  the  new  resistance  values  where  (2n  -   1)  is  the  number 
of  segments  in  the  curve. 

The  mechanization  of  the  switching  curve  involved  generation  of  the  functions, 

1  I     -  .       t  for  the  application  to  the  servo,  and  the  equivalent, 

m        ' 
yM         I    £     I       6     >  for  the  simulation.    In  the  first  case   6  was  multiplied  by 
2TmluJnl      u;n  ^ 

-rf  j/zrm      and  in  the  second  case    €  /u>   was  multiplied  by   ^/jl/zt     '. 

The  multiplication  conveniently  scaled  the  input,  e^,  and  also  reduced  the  problem 

to  one  of  generating    |ej  e±  for  either  case.    The  curve  was  approximated  by  five 

segments  from  -100  <e.  <  100  volts.    This  range  corresponds  to  that  from  minus 

one  to  plus  one  in  computer  units.    More  segments  were  tried,  but  the  change  in 

performance  of  the  servo  was  not  noticeable. 

The  steepest  slope  was  generated  by  one  amplifier  and  the  output  of  the  other 

amplifier  subtracted  from  it  to  generate  the  desired  curve.    The  slopes  used, 

range  of  input  voltages,  resistance,  and  bias  voltages  are  tabulated  below  for 

positive  e^  .    The  negative  half  of  the  curve  is  symmetric  with  respect  to  the 

origin.    A  one  megohm  feedback  resistor  is  used. 
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TABLE  H-l 
Diode  Function  Parameters  (No  Friction) 


ej(volts) 

Slope 

j  Change  of  Slope 

Rl    -    R2 

Bias  (Volts) 

0-10 

0 

i      ° 

" 

_ 

10-    40 

0.5 

0.5 

1.0 

5 

40  -  100 

1.25 

0.75 

|         0.67 

20 

! 

The  resulting  circuit  appears  in  Figure  n-3. 


-f  20 *   -20 "f 


FIGURE  H-3 
Diode  Function  Generator  (No  Friction) 
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The  outputs  of  the  two  amplifiers  are  shown  below  in  Figure  n-4  along 
with  the  contribution  of  the  -e .  input  to  amplifier  B. 
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FIGURE  H-4 
Diode  Function  (No  Friction) 


Although  this  particular  approximation  appears  very  crude,  a  comparison 
of  it  to  the  true  function  shows  there  is  no  more  than  one  volt  error  over  half 
the  range  and  four  volts  over  90%  of  the  range. 

By  adding  one  more  segment  on  each  end,  the  error  can  be  kept  within  one 
volt  over  97.  5%  of  full  scale. 


APPENDIX  IE 
MAGNETIC  AMPLIFIER  CHARACTERISTICS  AND  COMPENSATION 

To  determine  the  static  characteristics  of  the  magnetic  amplifier,  a 
control  current  was  supplied  from  an  operational  amplifier  and  the  output 
terminals  were  connected  to  a  1320  ohm  resistance  load.    This  load  represented 
ten  watts  at  rated  voltage.    The  static  response  curve  is  shown  in  Figure  ni-2. 

The  frequency  response  of  the  magnetic  amplifier  was  obtained  by  again 
driving  it  from  an  operational  amplifier  with  current  feedback.    The  operational 
amplifier  offered  no  bandwidth  restrictions  over  the  range  of  frequencies  con- 
cerned.   The  frequency  response  curve  in  Figure  ni-3  shows  that  the  output  was 
down  3  db  at  about  17  cycles  per  second.    By  adding  a  1.  6  jxi  condenser  and 
500  ohm  resistor  in  parallel  with  the  10k  resistor  to  ground,  the  output  was 
kept  quite  flat  out  to  about  200  cycles  per  second  after  which  it  falls  off  again. 
The  control  circuits  are  shown  in  Figure  D3-1. 
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FIGURE  m-1 
Magnetic  Amplifier  Control  Circuits 
(a)   Uncompensated   (b)   Compensated 
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